Probing Planck scale physics with IceCube 
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Neutrino oscillations can be affected by decoherence induced e.g. by Pianck scale suppressed inter- 
actions with the space-time foam predicted in some approaches to quantum gravity. We study the 
prospects for observing such effects at IceCube, using the likely flux of TeV antineutrinos from the 
Cygnus spiral arm. We formulate the statistical analysis for evaluating the sensitivity to quantum 
decoherence in the presence of the background from atmospheric neutrinos, as well as from plausible 
cosmic neutrino sources. We demonstrate that IceCube will improve the sensitivity to decoherence 
effects of Oi^E^ / Mp\) by 17 orders of magnitude over present limits and, moreover, that it can probe 
decoherence effects of 0{E'^ / Mpi) which are well beyond the reach of other experiments. 
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I. INTRODUCTION 

Despite many decades of intense effort, a satisfactory 
theory of quantum gravity is yet to see the light of day. 
Moreover, it is generally thought that the quantum ef- 
fects of gravity may never be experimentally accessible 
because they would be manifest only at the Planck scale, 

Afpi = ~ 1.2 X 10^^ GeV. However, gravity, be- 

ing a non-renormalizable interaction in the language of 
quantum field theory, may leave a distinctive imprint at 
energies much lower than the Planck scale if it violates 
some fundamental symmetry of the effective low energy 
theory, akin to the violation of parity in nuclear radioac- 
tive decay at energies far below the true scale of the re- 
sponsible weak interaction 0. For example, if quantum 
space-time has a 'foamy' structure in which Planck length 
size black holes form and evaporate on the Planck time 
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scale, then there may be a loss of quantum information 
across their event horizons, providing an 'environment' 
that can induce decoherence of apparently isolated mat- 
ter systems 

The particle most sensitive to such effects would ap- 
pear to be the neutrino because oscillations between neu- 
trino flavours are a pure quantum phenomenon in which 
the density matrix, p, has the properties of a projection 
operator: Tt — Tr p — 1. A heuristic view of decoher- 
ence induced by interactions with the virtual black holes 
in the space-time foam is as follows. Since black holes are 
believed not to conserve global quantum numbers, neu- 
trino flavor is randomized by interactions with the virtual 
black holes. The result of many interactions then is to 
equally populate all three possible flavors. Because black 
holes do conserve energy, angular momentum, and elec- 
tric and color charge (unbroken gauged quantum num- 
bers), the neutrino interacting with the virtual black 
hole does re-emerge as a neutrino. In this connection, 
it has been noted already that the results from the 
Super-Kamiokande atmospheric neutrino experiment \m 
and the K2K long baseline oscillation experiment 
interpreted in terms of a 2-generation <-> Uj- flavor 
transition, can probe decoherence effects with high sen- 
sitivity, supplementing laboratory tests based on K'^K 
oscillations and neutron interferometry P|. 
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It has recently been suggested that antineutrinos 
originating in the decay of neutrons from candidate cos- 
mic ray sources in the Galaxy can provide an even 
more sensitive probe. The effects of quantum decoher- 
ence would alter the flavor mixture to the ratio, v,, : v^^ : 
Vr — 1 : \ 1, regardless of the initial flavor content. 
Since decoherence effects are weighted by the distance 
travelled by the (anti) neutrinos, this means that if a V- 
flux with ratio of flavors ^ 1 : 1 : \ were to be observed 
from an astrophysical source, then strong constraints can 
be placed on the energy scale of quantum decoherence, 
surpassing current bounds by over 10 orders of magni- 
tude. However if a 1:1:1 ratio is observed, this will not 
imply that quantum decoherence is responsible, since the 
dominant source of the (anti) neutrinos may simply not 
be neutron /3-decay. In this paper, we pursue this idea 
further and formulate the statistical analysis necessary 
for obtaining bounds on quantum decoherence from ex- 
pected future detections of cosmic neutrinos. 

In Sec.^ we identify a candidate neutron source in the 
vicinity of the Earth: Cygnus 0B2. We review all exist- 
ing data on the Cygnus region and show that observed 
directional signals at high energies 0, 0| can plausibly 
be ascribed to a neutron source with an energy spectrum 
oc E:^^-^. In particular, because of neutron decay, the 
expected anisotropy is well below limits reported by the 
CASA-MIA ^] and KASCADE [ij experiments. We 
summarize the estimate of the corresponding antineu- 
trino flux In Sec. IIIII we discuss the effects of de- 
coherence on high energy neutrino propagation adopt- 
ing the quantum dynamical semi-group formalism, where 
the Lindblad operators |l3j | describe (anti) neutrino cou- 
plings to the space-time foam. In Sec. IIVI we estimate 
in detail the sensitivity of the IceCube detector |3| to 
the F-Cygnus beam and its ability to constrain the ef- 
fects of quantum decoherence. We evaluate the signal- 
to-noise for both track and shower events, taking into 
account the atmospheric neutrino background, as well 
as a possible contribution to the neutrino flux from the 
TeV 7-ray source TeV J2032-f4130, recently discovered 
by the HEGRA experiment jl^ il^] in the direction of 
the Cygnus spiral arm. Armed with these event rates, 
we formulate the statistical analysis required to study 
the sensitivity to quantum decoherence effects. We show 
that IceCube can improve the sensitivity over present 
probes of decoherence by 4 to 17 orders of magnitude, and 
moreover, that it is sensitive to strongly energy depen- 
dent decoherence effects suppressed by multiple powers 
of the Planck scale which are beyond the reach of other 
experiments. Finally, in Sec. we confront our results 
with theoretical suggestions for quantum decoherence. 



II. 



ANTINEUTRINOS FROM CYGNUS OB2 



Massive star forming regions are the engines of star- 
burst galaxies. They generate large numbers of UV pho- 
tons which ionize the interstellar medium, as seen in 
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Cygnus 0B2 Field: HEGRA CT-System 
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FIG. 1: Skymap of correlated event excess significance (a) 
from all HEGRA lACT-System data (3° x 3° FoV) centered 
on the TeV source J2032+4130. Nearby objects are also 
shown: 95% contours for 3 EGRET sources (indicated by 
the ovals), their possible X-ray associated counterparts (as 
given in Ref. [l^), and Cygnus X-3. The center of gravity 
with statistical errors and intrinsic size (standard deviation 
of a 2-dim Gaussian, asrc) are indicated by the white cross 
and white circle, respectively. The TeV source, J2032-(-4130, 
is positioned at the edge of the error circle of the EGRET 
source 3EG J2033-I-4118, and within the core circle of the 
extremely dense OB stellar association Cygnus OB2 Hi. 



microwave, radio, and Ha recombination line emission. 
They are important sources of interstellar dust heating 
which results in significant infrared emission. The strong 
winds of their massive O stars, which should pass through 
the Wolf-Rayet phase and explode as supernovae, release 
considerable amounts of kinetic energy creating a rarefied 
hot (~ 10^ K) superbubble that emits X-rays; its cavity 
can eventually be discerned from observation of H i and 
CO lines of the interstellar gas. 

Such regions are also likely sites for cosmic ray acceler- 
ation. The massive stars synthesize considerable amounts 
of heavy nuclei that are released either by stellar winds or 
during the subsequent supernova explosions. Moreover, 
the young stellar population can create time-correlated, 
clustered supernova remnants, where through coopera- 
tive acceleration processes the energy of the accelerated 
nuclei can be boosted above the '^10^ GeV cutoff of indi- 
vidual remnants. Thus the usual Fermi mechanism might 
be able to accelerate cosmic rays all the way up to the 
'ankle', where the steeply falling (oc E^^'^) cosmic ray 
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spectrum flattens to E^^^ 0|- An immediate conse- 
quence of this nucleus-dominance picture is the creation 
of free neutrons through the photodisintegration of the 
nuclei on the intense ambient photon fields. 

Independent evidence may be emerging for such a cos- 
mic ray accelerator in the massive star forming region 
Cygnus 0B2, a cluster of several thousands of hot young 
OB stars with a total mass of 10^ Mq [l^- At a rel- 
atively small distance to Earth {d « 1.7 kpc), this is 
the largest known stellar association, with a diameter 
of ~ 60 pc and a core radius of ~10 pc. The cluster 
age has been estimated from isochrone fitting to be 



3-4 Myr, where the age range may reflect a non-coeval 
star forming event. The HEGRA experiment has de- 
tected an extended TeV 7-ray source, J2032+4130, on 
the outer edge of Cygnus 0B2 with no clear counterpart 
and a spectrum which can be modelled in terms of either 
hadronic or leptonic processes However the fail- 

ure of Chandra and VLA to detect significant levels of 
X-rays or radiowaves which would signal the presence of 
high energy electrons .20] argues for a hadronic mecha- 
nism. Above 1 TeV, the HEGRA data can be fitted by a 
simple power law 
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The model proposed |2l| is that protons are accelerated 
by the collective effects of stellar winds from massive O 
and B stars and only the high energy particles penetrate 
and interact in the innermost dense parts of the winds. 
The colliding protons generate tt^'s which produce the 
observed 7-rays. Convection prevents low energy protons 
from entering the dense wind region thus explaining the 
absence of MeV-GeV photons. 

At very high energies {E > 10^'^ GeV) evidence has 
also been presented for neutral particles from the Cygnus 
spiral arm. AGASA has found a 4a correlation of the 
arrival direction of cosmic rays at these energies with 
the Galactic Plane (CP) jH. The CP excess, which is 
roughly 4% of the diffuse flux, is mostly concentrated in 
the direction of the Cygnus region Evidence at the 
3.2cr level for a CP enhancement at similar energies has 
also been reported by the HiRes Collaboration |2J| . The 
primary particles must be neutral (and stable) in order to 
preserve direction while propagating through the galactic 
magnetic field. In principle they can be photons but this 
is hard to reconcile with the complete isotropy observed 
up to - 10^-^ GeV by KASCADE Intriguingly, 
time-dilated neutrons can reach the Earth from typical 
Galactic distances when their energy exceeds ~10^ GeV 
so it is reasonable to ask whether these might in fact be 
the primaries. 

The GP anisotropy is observed over the energy range 
10^-^ to lO^ '^ GeV. The lower threshold requires that 
only neutrons with energy > 10^ GeV have a boosted 
CT„ sufficiently large to serve as Galactic messengers. 
The decay mean free path of a neutron is c r„ r„ = 
9.15 {En/lO^ GeV) kpc, the hfetime being boosted from 
its rest-frame value, r„ — 886 seconds, to its lab value by 
Pn = En/rrin- Actually, the broad scale anisotropy from 
the direction of the GP reported by Fly's Eye j2J| peaks 
in the energy bin 10^'^ — 10^ GeV, but persists with less 
significance to energies as low as 10^'^ GeV. This implies 



that if neutrons are the carriers of the anisotropy, there 
needs to be some contribution from at least one source 
closer than ^2 kpc. Interestingly, the full Fly's Eye 
data includes a directional signal from the Cygnus re gion 
which was somewhat lost in unsuccessful attempts ^ [l^l 
to relate it to 7-ray emission from Cygnus X-3. As shown 
in Fig. n Cygnus 0B2 is very close to the line of sight to 
Cygnus X-3, which is in fact ~ 8 kpc farther away than 
the stellar association. 

The upper cutoff reflects an important feature of pho- 
todisintegration at the source: heavy nuclei with ener- 
gies in the vicinity of the ankle will fragment to neutrons 
with energies about an order of magnitude smaller. To 
account for the largest neutron energies, it is necessary 
to continue the heavy nucleus spectrum to energies above 
the ankle Note that the emerging harder extragalac- 
tic spectrum will overwhelm the steeply falling galactic 
population at these energies. In order to fit the spectrum 
in the anisotropy region and maintain continuity to the 
ankle region without introducing a cutoff, the AGASA 
Collaboration required a spectrum cx E^^ or steeper ii2^ . 

In what follows, we model the neutron spectrum with 
a single power law refiecting the average shape of the 
diffuse cosmic ray spectrum between 10^ and 10*'^ GeV, 
specifically: 
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By integrating the spectrum between Ei — 10*'^ GeV 
and E2 = 10^'^ GeV l2a|, we can normalize to the ob- 
served integrated fiux |2^ : 
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which yields C = 1.15 x 10^" km ^ yr ^. We emphasize 
again that the neutron primaries hypothesis predicts a 
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significant signal above the diffuse cosmic ray flux only 
at energies > 10^-^ GeV. Figure |21 shows the damping 
due to neutron decay which attenuates the directional 
signal at low energies. The predicted damped signal for 
a source at 1.7 kpc is well below direct limits from the 
CASA-MIA [m and KASCADE [l^ experiments. 

For every surviving neutron at ~ 10^ GeV, there are 
many neutrons at lower energies that do decay via n — > 



p + e~ + Ve. The proton is bent by the Galactic magnetic 
field and the electron quickly loses energy via synchrotron 
radiation, but the Ve travels along the initial neutron 
direction, producing a directed TeV energy beam which 
is potentially observable. 

The basic formula that relates the neutron flux at the 
source to the antineutrino flux observed at Earth is M: 
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The variables appearing in Eq. are the antineutrino 
and neutron energies in the lab (£77, i?„), the antineu- 
trino angle with respect to the direction of the neutron 
momentum in the neutron rest-frame {6^7), and the an- 
tineutrino energy in the neutron rest-frame (eu). The 
last three variables are not observed by a laboratory 
neutrino-detector, and so are integrated over. The ob- 
servable £"77 is held fixed. The delta-function relates the 
neutrino energy in the lab to the three integration vari- 
ables. The parameters appearing in Eq. Q are the neu- 
tron mass and rest-frame lifetime (mn and Tn)- Finally, 
dP/dt-jj is the normalized probability that the decaying 
neutron produces an antineutrino with energy exr in the 
neutron rest-frame. Note that the maximum antineu- 
trino energy in the neutron rest frame is very nearly 
Q = jTin - nip - nie = 0.71 MeV. 

The integral neutrino flux, F-^{> Ejj) = /g_d£V 
is particularly useful for experiments having a neutrino 
detection-efficiency that is independent of neutrino en- 
ergy, or nearly so. Our calculated integral flux, normal- 
ized to the integrated neutron flux in Eq. (O , is shown in 
Fig. 13 Note that the nuclear photodisintegration thresh- 
old implies an infrared cutoff on the primary neutron 
energy at the source, which in turn leads to a low energy 
cutoff of O(TeV) on the integral flux. 



III. DECOHERENCE EFFECTS IN HIGH 
ENERGY NEUTRINO PROPAGATION 
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FIG. 2: The integrated neutron flux expected from Cygnus 
OB2 (dashed-dotted line) is superimposed over the integrated 
fluxes observed from the Cygnus region by the Fly's Eye |3 
and AGASA ITol . experiments. Also shown is the 7-ray 
flux reported by the HEGRA experiment and the upper 
limits on neutral particles derived from the CASA-MIA pdl| 
and KASCADE dj experiments. The solid line is a fit to 
the HEGRA data and the dashed line is the extrapolation to 
unobserved energies. 



Even though the flux of antineutrinos produced by 
Cygnus 0B2 is pure Fg, the antineutrinos observed at 
Earth will be distributed over all flavors. This is because 
of neutrino oscillations, as well as possible decoherence 
effects induced over long distances (more on this below) . 
In the standard treatment of neutrino oscillations, neu- 
trino flavor eigenstates, a — e,ii,T, are expanded 

in terms of mass eigenstates, \ vi), i — 1,2,3, through 

3 

a (unitary) matrix, U, defined by {va) = Vj). 



This implies that the density matrix of a flavor state, 
p", can be expressed in terms of mass eigenstates by 
= Wa)Wa\ = J2t,jU*.^Uaj\iyi){iyj\. This Is a pure 
quantum system, therefore the density matrix satisfies 
Tr = Tr p = 1. To get the transition amplitude, we 
evolve the system quantum mechanically with the Liou- 
ville equation 



dp 
dt 



-I [H, p] 



(5) 



5 




E [GeV] 

FIG. 3: Integrated flux of 77^ + 77e + v-r (solid line) predicted 
to arrive at Earth from the direction of the Cygnus region. 
Also shown are the integrated Vfj, + and Ve + Ve atmo- 
spheric fluxes for an angular bins of 1° x 1° and 10° x 10°, 
respectively. The shaded band indicates the region excluded 
by the AMANDA experiment I^^.The fluxes of neutrinos in- 
ferred from HEGRA measurements of the 7-ray flux are also 
shown: the lower line is based on the assumption of p7 in- 
teractions, whereas the upper line is based on pp interactions 
(the charged/neutral pion-production ratio depends on the in- 
teraction). In each case the solid portion of the line indicates 
the region where HEGRA data is available and the dashed 
part is an extrapolation to unobserved energies. 



where H is the Hamiltonian of the system. For 
5rnljd/2E ^ 1, the phases induced by the mass spht- 
ting, drriij = ml — m| , will be erased by uncertainties in 
d and E, yielding for the transition probability between 
flavor states a and (3 [29l |: 

= S^0-2j2 Re Up, U^j (7* J . (6) 

The prediction for the flavor population at Earth due to 
standard flavor-mixing (i.e. with no spacetime dynam- 
ics) of the pure Ve source is J2j I Uej \ ^ | Uaj P ~ 5 1 C/q2 | ^ + 
||C^aiP for flavor a, which leads to the flavor ratios 
^5:2:2. This is very different from the democratic 
1:1:1. 

The Hamiltonian evolution in Eq. |SJ is a character- 
istic of physical systems isolated from their surround- 
ings. The time evolution of such a quantum system is 
given by the continuous group of unitarity transforma- 
tions, Ut — e~'^*, where t is the time. The existence of 
the inverse of the infinitesimal generator, which is a 
consequence of the algebraic structure of the group, guar- 
antees reversibility of the processes. For open quantum- 



mechanical systems, the introduction of dissipative ef- 
fects lead to modifications of Eq. jSJ that account for the 
irreversible nature of the evolution. The transformations 
responsible for the time evolution of these systems are 
defined by the operators of the Lindblad quantum dy- 
namical semi-groups 13]. Since this does not admit an 
inverse, such a family of transformations has the property 
of being only forward in time. 

The Lindblad approach to decoherence does not re- 
quire any detailed knowledge of the environment. The 
corresponding time evolution equation for the density 
matrix takes the form: 

^ = -i[H,f,,p]+V[p] , (7) 

where the decoherence term is given by 

^W = -^E([^J'^^^] + [&.P,6]]) • (8) 
j 

Here, Hcs = H + is the effective Hamiltonian of the 
system, H is its free Hamiltonian, accounts for possi- 
ble additional dissipative contributions that can be put in 
the Hamiltonian form, and {6^ } is a sequence of bounded 
operators acting on the Hilbert space of the open quan- 
tum system, and satisfying b'^-bj e B{J-L), where 
B{T-L) indicates the space of bounded operators acting on 

n. 

The intrinsic coupling of a microscopic system to the 
space-time foam can then be interpreted as the exis- 
tence of an arrow of time which in turn makes possi- 
ble the connection with thermodynamics via an entropy. 
The monotonic increase of the von Neumann entropy, 
S{p) = — Tr (p Inp), implies the hermiticity of the Lind- 
blad operators, bj = 5^ jSQi]. In addition, the conservation 
of the average value of the energy can be enforced by tak- 
ing [H, bj] =0 [33. The Lindblad operators of an TV-level 
quantum mechanical system can be expanded in a basis 
of matrices satisfying standard commutation relations of 
Lie groups. For a 3-level system, the basis comprises the 
eight Gell-Mann SU(3) matrices plus the 3x3 identity 
matrix [32i |. 

As mentioned above, the theoretical approach provided 
by Lindblad quantum dynamical semi-groups is a very 
general in the sense that no explicit hypothesis has been 
made about the actual interactions causing the loss of 
coherence. Following Ref. [s^, we adopt an expansion 
in a 3 flavor basis with a diagonal form for the 9x9 
decoherence matrix, V. Note that neutrinos oscillate 
among flavors separately between particle and antiparti- 
cle sectors and so the respective decoherence parameters 
for antineutrinos can be different from the correspond- 
ing ones in the neutrino sector. Upon averaging over the 
rapid oscillation for propagation between Cygnus 0B2 
and the Earth, only the diagonal Gell-Mann matrices sur- 
vive, and so the transition probability for antineutrinos 
takes the form [ssf : 
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where 73 and 7g are eigenvalues of the decoherence ma- 
trix for antineutrinos. Note that in Eq. we set the 
CP violating phase to zero, so that all mixing matrix 
elements become real. Furthermore, under the assump- 
tions that CPT is conserved and that decoherence effects 
are negligible at present experiments, the values of the 
mixing angle combinations appearing in Eq. ||5J) can be 
well determined by the usual oscillation analysis of solar, 
atmospheric, LBL and reactor data 0|. In what follows, 
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we will assume that CPT is conserved both by neutrino 
masses and mixing as well as in decoherence effects. Note 
however that since the decoherence effects in the present 
study primarily affect antineutrinos, our results will hold 
for the antineutrino decoherence parameters exclusively 
if CPT is violated only through quantum-gravity effects. 

Now, we require further 73 = = 7 ( = 73 = under 
CPT conservation) so that Eq. can be rewritten for 
the case of interest as: 
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We make this simplification only to emphasize the pri- 
mary signature of quantum decoherence, namely that af- 
ter travelling a sufhciently long distance the flavor mix- 
ture is altered to the ratio 1:1:1, regardless of the 
initial flavor content. Consequently, if a flux of antineu- 
trinos were to be observed from the Cygnus spiral arm 
with a flavor ratio 7^ 1 : 1 : 1, strong constraints can be 
placed on the decoherence parameter 7. 

Using the results of the up-to-date oscillation anal- 
ysis of solar, atmospheric, LBL and reactor data j3J] we 
obtain the following values and 95% confidence ranges 



U.iUD_o.o82 

n 100-1-0.089 
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The numbers given in Eq. are obtained using the 
same techniques as described in Ref. ^34] but including 
the final SNO salt phase data [ssf . 

We assume a phenomenological parametrization for 
the eigenvalues of the decoherence matrix 



7 = /«„ (£;,/GeV)", 



(12) 



with the integer n 6 [—1, 3]. This allows a straightfor- 
ward comparison with existing limits. Equipped with 
Eqs. ((Tn|l . (fTTT) and ((T^ . we now proceed to determine 
the sensitivity of IceCube to decoherence effects. 



IV. SENSITIVITY REACH AT ICECUBE 

In deep ice, neutrinos are detected through the obser- 
vation of Cerenkov light emitted by charged particles pro- 
duced in neutrino interactions. In the case of an incident 
high-energy muon neutrino, for instance, the neutrino in- 
teracts with a hydrogen or oxygen nucleus in the deep ice 
and produces a muon traveling in nearly the same direc- 
tion as the neutrino. The blue Cerenkov light emitted 
along the muon's kilometer-long trajectory is detected 
by strings of PhotoMultiplier Tubes (PMTs) deployed 
at depth shielded from radiation. The orientation of the 
Cerenkov cone reveals the neutrino direction. There may 
also be a visible hadronic shower if the neutrino is of suf- 
ficient energy. 

The Antarctic Muon And Neutrino Detector Array 
(AMANDA) [35 , using natural 1 mile deep Antarctic ice 
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as a Cerenkov detector, has operated for more than three 
years in its final configuration: 19 strings instrumented 
with 680 PMTs. IceCube the successor experiment 
to AMANDA, is now under construction. It will consist 
of 80 kilometer-length strings, each instrumented with 60 
PMTs spaced by 17 m. The deepest module is 2.4 km be- 
low the surface. The strings are arranged at the apexes of 
equilateral triangles 125 m on a side. The instrumented 
(not effective!) detector volume is a full cubic kilome- 
ter. A surface air shower detector, IceTop, consisting 
of 160 Auger-style 0| Cerenkov detectors deployed over 
1 km^ above IceCube, augments the deep-ice component 
by providing a tool for calibration, background rejection 
and air-shower physics. The expected energy resolution 
is ±0.1 on a logio scale. Construction of the detector 



started in the Austral summer of 2004/2005 and will con- 
tinue for 6 years, possibly less. At the time of writing, 
data collection by the first string has begun. 

At IceCube, the events are grouped as either muon 
tracks or showers. Tracks include muons resulting from 
both cosmic muons and from Charged Current (CC) in- 
teraction of muon neutrinos. The angular resolution 
for muon tracks « 0.7° 38] allows a search window of 
1° X 1°. This corresponds to a search bin solid angle of 
Afiioxio « 3 X 10~^ sr. In order to reduce the back- 
ground from cosmic muons, we adopt here the quality 
cuts referred to as "level 2" cuts in Ref. [T^ . 

In our semianalytical calculation |39|. we estimate the 
expected number of induced tracks from the Cygnus 
0B2 source antineutrino flux as 
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where 



(14) 



is the differential antineutrino flux which arrives at the 
Earth with flavour a after oscillation of the Ve in Eq.Q). 



Here 



do-cc 
dEO 



{E-ij,E'^) is the differential CC interaction 

cross-section producing a muon of energy E'^, rtx is the 
number density of nucleons in the matter surrounding 
the detector, and T is the exposure time of the detector. 
After being produced, the muon propagates through the 
rock and ice surrounding the detector and loses energy. 
We denote by F{E^, E^'^,1) the function that represents 
the probability of a muon produced with energy E'JJ, ar- 
riving at the detector with energy E^'^, after traveling a 
distance, I. The details of the detector are encoded in 
the effective area, We use the parametrization of 

the A^g described in Ref. ^3 to simulate the response 
of the IceCube detector after events that are not neutri- 
nos have been rejected (this is achieved by quality cuts 
referred to as "level 2" cuts in Ref. 111). The" mini- 
mum track length cut is Zmin ~ 300 m and we account 
for events with E^" > 500 GeV. With this we find that, 
assuming standard neutrino oscillations, one expects a 



total of 212 X Pjy^^-p^ — 48 T7^-induced tracks from the 
Cygnus 0B2 source in 15 years of observation. 

Showers are generated by neutrino collisions — 
i^e or CC interactions, and all Neutral Current (NC) 
interactions — inside of or nearby the detector, and by 
muon bremsstrahlung radiation near the detector. For 
showers, the angular resolution is significantly worse than 
for muon tracks. In our analysis, we consider a shower 
search bin solid angle, ASliooxio"- Normally, a reduction 
of the muon bremsstrahlung background is effected by 
placing a cut of 40 TeV on the minimum reconstructed 
energy For Cygnus 0B2, this strong energy cut 

is not needed since this background is filtered by the 
Earth. Thus we account for all events with shower energy 
Esh > = 1 TeV. The directionality requirement, 

however, implies that the effective volume for detection 
of showers is reduced to the instrumented volume of the 
detector, 1 km^, because of the small size of the showers 
(less than 200 m in radius) in this energy range. 

We can now calculate the expected number of showers 
from the Cygnus 0B2 source to be: 



■sh,NC 



(15) 



where 



■sh,CC 



(16) 



and 



iL rsh.NC T-T -» 1 



E--E 



•.h -^sh a=e,fj,,T 



(17) 
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Here, ^^{Etj,EL) is the differential NC interaction 
cross section producing a secondary antineutrino of en- 
ergy, E^. In writing Eqs. H16(l and (|17|l we are assum- 
ing that for contained events the shower energy corre- 
sponds with the interacting or Vj- antineutrino energy 
{Esh = E^r) in a CC interaction, while for NC the shower 
energy corresponds to the energy in the hadronic shower 
Esh = -E57 — EL = E-jy-y, where y is the usual inelasticity 
parameter in DIS. In total, within the framework of stan- 
dard oscillations, we expect 25 showers from the Cygnus 
0B2 source in 15 years of operation. 

We now turn to the estimate of the background. There 
are two different contributions — atmospheric neutrinos 
and additional fluxes of extraterrestrial neutrinos. For 
the "conventional" atmospheric neutrino fluxes arising 
from pion and kaon decays, we adopt the 3-dimensional 
scheme estimates of Ref. |41|, which we extrapolate to 
match at higher energies the 1-dimensional calculations 
by Volkova 42|. We also incorporate "prompt" neutri- 
nos from charm decay as calculated in Ref. We 
obtain the number of expected track and shower events 
from atmospheric neutrinos as in Eqs (|13|l , (|16|l , and H17|) 

with — (^i/) being the i^e and i/^ atmospheric neu- 
trino fluxes integrated over a solid angle of of 1° x 1° (for 
tracks) and 10° x 10° (for showers) width around the di- 
rection of the Cygnus 0B2 source {0 = 131.2°). We get 
an expected background of 14 atmospheric tracks and 47 
atmospheric showers in 15 years. Of the 47 showers, 16 
correspond to CC interactions while 31 correspond to 
i^fj, NC interactions. The large yield of NC events is due 
to the fact that at these energies, the atmospheric flux 
contains a very unequal mix of neutrino flavors (with ra- 
tios « 1 : 20 : 0). We have also verified that if we increase 
the minimum shower energy cut to 5 TeV, I'e CC and i/^ 
NC contribute in equal amounts to the number of atmo- 
spheric showers. This is in agreement with simulations 
by the AMANDA Collaboration 44]. 

We turn now to the discussion of background events 
from other extraterrestrial sources. As discussed in 
Sccini the TeV 7-ray flux reported by the HEGRA Col- 
laboration [1^ in the vicinity of Cygnus 0B2 is likely to 
be due to hadronic processes: the 7-rays can be directly 
traced to the decay tt^'s produced through inelastic pp 
collisions Since 7r°'s, tt+'s, and 7r~'s are produced 

in equal numbers, we expect two photons, two i^e's, and 
four Vf^'s per tt". On average, the photons carry one- 
half of the energy of the pion, and the neutrinos carry 
one-quarter. During propagation, the i/^'s will partition 
themselves equally between i/^'s and i^t's on lengths large 
compared to the oscillation length and so one finds at 
Earth a nearly identical flux for the three neutrino fla- 
vors I45II: 



c\F dF 

''^J^^E,^E^I2)^2^{E^ 



dE^ 



dE^ 



(18) 



For p7 — > N-K interactions, it can easily be shown us- 
ing the A-approximation that the resulting neutrino flux 



is about a factor of 4 smaller |4g . For the purposes of 
setting an upper bound on the neutrino flux we ignore 
all other sources near J2032-I-4130 because their steady 
emission in 7-rays is estimated to be smaller by more over 
a factor of 5 than the source of interest |^ . Substituting 
Eq. into Eq. H18|l we obtain the corresponding back- 
ground from neutrinos with flavor ratios 1:1:1. As can 
be seen in Fig.|3| the background is dominated by atmo- 
spheric neutrinos. Thus after 15 years of data collection 
we expect 18 tracks and 1 shower from J2032-I-4130 for 
standard oscillations. In summary, the directional beam 
from the Cygnus region provides a statistically significant 
signal to probe anomalous oscillations in the antineutrino 
sector. 

We will now discuss how to isolate the possible sig- 
nal due to decoherence in the antineutrinos from Cygnus 
0B2 from the atmospheric background and possible fluc- 
tuations in the event rate due to unknown diffuse fluxes 
of extraterrestrial neutrinos. In general, we can predict 
that the expected number of track and shower events in 
the direction of the Cygnus 0B2 source to be 



-MY 



ATM 



SS ; 



(19) 
(20) 



The first term corresponds to antineutrinos from neu- 
tron /9-decay. In the presence of decoherence effects these 
event rates can be computed from Eqs. l(T^ . ((TB|l and l(T7|) 
with flavour transition probabilities given in Eq. IjlOl) 
with c? = 1.7 kpc. The second term refers to atmospheric 
(anti) neutrinos (A/Xtm = 14, A/^tm — 47 for 15 years 
of exposure). The third term takes into account addi- 
tional contributions from a diffuse flux of (anti)ncutrinos 
produced via charged pion decay. In principle, decoher- 
ence effects may also affect the expected number of events 
from this diffuse flux. However given that the flavour ra- 
tios both from oscillation and decoherence are very close 
to 1:1:1 for the case of neutrinos produced via charge 
pion decay, we find that there is no difference in the sen- 
sitivity region if decoherence effects are included or not 
in the evaluation of Agg and A/^g. They are normalized 
to the maximum expected flux from J2032+4130 by a 
factor X = 7V^yi8 = A^||/l. 

Altogether, the quantities A/"'"^ and M^^, as defined in 
Eqs. H19|) and H20|) . can be regarded as the theoretical 
expectations of these events rates, corresponding to dif- 
ferent points in the x — Kn parameter space. For a given 
set of observed rates, ^^'^ -^obs' ^'^^ curves are ob- 
tained in the two-dimensional parameter space by setting 
Af^s = and A/'^Jj^, = A/"*"^. These curves intersect at a 
point, yielding the most probable values for the flux and 
decoherence scale for the given observations. Fluctua- 
tions about this point define contours of constant in 
an approximation to a multi-Poisson likelihood analysis. 
The contours are defined by ,4^ 



sh, tr 



AA^-A-^bs+ACb^ In 



obs 



AO 
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levels currently probed, we now comment briefly on the 
theoretical implications. 

Any type of high energy/short distance space-time 
foam interaction given in Eq. can be understood in 
analogy with the tracing out of the degrees of freedom 
of a thermal bath (with temperature T) with which the 
open system (in our case a neutrino beam) interacts. A 
simple version of the interaction with the bath can be 
written as ifint = 'Ylij ^ji'^ + 0.^)1 where a, at are raising 
and lowering operators for space-time foam excitations, 
with (a^a) = (e^*--"-/^ - 

The energy behavior of 7 depends on the dimensional- 
ity of the operators bj. But care must be taken, since V is 
bilinear in the bj, and due to the hermiticity requirement, 
each bj is itself at least bilinear in the neutrino fields ip. 
Examples are 



(27) 



bj cx / d'^x (idtYii 



FIG. 4: IceCube's sensitivity to quantum decoherence assum- 
ing 7V*bg = 62 tracking events and N^^^ — 72 showering events 
from Cygnus. The regions above and to the right of the solid, 
dashed and dotted lines can potentially be excluded at 90%, 
95%, and 99% confidence level, respectively. 



As illustration, in Fig. 0] we show for the case n = 0, 
the expected constraints on at 90, 95 and 99% CL 
for 2 d.o.f. if observations turn out to be in agreement 
with standard neutrino oscillation expectations, taking 
Mll^ = 62 and M^^^ = 72 (and no diffuse flux). Similar 
regions can be obtained for other choices of n. 

Marginalizing with respect to x, we extract the fol- 
lowing 1 degree-of-freedom bounds on the decoherence 
parameters 
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at 90 (99) % CL. These should be compared with the 
current 90% CL upper limits on the decoherence pa- 
rameter from the Super-Kamiokande and K2K data: 
K_i < 2.0 X 10-21 GeV , Ko < 3.5 x lO^^s QeV and 
K2 < 9.0 X 10-28 GeV y. It is clear that IceCube will 
provide a major improvement in sensitivity to the possi- 
ble effects of quantum gravity. 



DISCUSSION 



where j = 0,l,2....A Fourier expansion of the fields 
'0, V'^j inserted into Eq. (|27|l . gives the energy behavior 
bj cx iii, and hence 7 cx E^^ . This restriction of the en- 
ergy behavior to non-negative even powers of Ej, may 
possibly be relaxed when the dissipative term is directly 
calculated in the most general space-time foam back- 
ground. 

An interesting example is the case where the dis- 
sipative term is dominated by the dimension-4 opera- 
tor 61, J X ip^ idtip , yielding the energy dependence 
7 cx E^/Mpi. This is characteristic of non-critical string 
theories where the space-time defects of the quantum 
gravitational "environment" are taken as recoiling D- 
branes, which generate a cellular structure in the space- 
time manifold 50] . 

Although the cubic energy dependence 7 cx is not 
obtainable from the simple operator analysis presented 
above, it may be heuristically supported by a general 
argument that each of the bj must be suppressed by at 
least one power of Mpi, giving a leading behavior 



j = k El/M, 



PI • 



(28) 



Here k is a dimensionless parameter which by natural- 
ness is expected to be 0(1). Decoherence effects with 
this energy behavior are undetectable by existing exper- 
iments. However, since the loss of quantum coherence is 
weighted by the distance travelled by the antineutrinos, 
by measuring the Iz-Cygnus beam IceCube will attain a 
sensitivity down to k < 3.0 x 10"^ at 99% CL, well below 
the natural expectation. 

Finally, we note an interesting aspect of the k_i limit. 
For n = —1, a non- vanishing 7 in Eq. I|10l) can be related 
to a finite Ve lifetime in the lab system 51|: 



-d/ri. 



-dniTj^/ETj 



(29) 



Having demonstrated that IceCube will be able to set 
bounds on quantum decoherence effects well beyond the 



where r^y^ is the antineutrino rest frame lifetime and rriTi^ 
its mass. Therefore the 90% bound from IceCube on k_i 
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can be translated into 

^ > lO^-* GeV-^ = 6.5 s eV"^ . (30) 

This corresponds to an improvement of about 4 orders of 
magnitude over the best existing bounds from solar neu- 
trino data jHll , and of course gives results comparable to 
the reach derived previously for neutrinos decaying over 
cosmic distances jS^- It should be noted that although 
the similar algebraic structure of the decoherence term in 
Eq. (|10|l and a decaying component in the neutrino beam 
provide a bound on the neutrino lifetime, these are con- 
ceptually two different processes. The decoherence case 
can be viewed as a successive absorption and re-emission 
of a neutrino from the quantum foam with change in fla- 
vor but no change in the average energy because of the 
condition [H, bj] = 0. This contrasts with the decay pro- 
cess which involves the emission of an additional particle. 

In conclusion, the IceCube experiment will be sensi- 
tive to the effects of quantum decoherence at a level well 
beyond current limits. In particular, antincutrinos pro- 
duced in the decays of neutrons from Cygnus 0B2 pro- 
vide an excellent source in which to search for such ef- 
fects [E^]. Although the precise conclusions depend on 
the model considered, we find that in general IceCube 



can achieve a sensitivity of more than 10 orders of mag- 
nitude beyond current bounds on decoherence through 
observations of Galactic sources of neutrons. 
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